Global culture: A noise induced transition in finite systems 
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We analyze the effect of cultural drift, modeled as noise, in Axelrod's model for the dissemination of culture. 
The disordered multicultural configurations are found to be metastable. This general result is proven rigorously 
in d = 1, where the dynamics is described in terms of a Lyapunov potential. In d = 2, the dynamics is 
governed by the average relaxation time T of perturbations. Noise at a rate r < T _1 induces monocultural 
configurations, whereas r > T _1 sustains disorder. In the thermodynamic limit, the relaxation time diverges 
and global polarization persists in spite of a dynamics of local convergence. 
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Concepts and tools developed in the framework of Statis- 
tical and Nonlinear Physics have been shown to be useful in 
identifying general mechanisms behind collective behavior in 
social dynamics [[!]]. Beyond the use of equilibrium ideas 
based on minimization of a potential or optimization meth- 
ods [[|], models of interacting social agents, as for example 
generalized voter models of opinion formation [|J, fall within 
the context of the study of nonequilibrium phase transitions in 
lattice models [B]. More generally, these are stochastic spatial 
models also considered in population dynamics or evolution- 
ary biology 

Rather independently of this framework, but closely related 
in spirit, is the model put forward by Axelrod to describe the 
dissemination of culture among interacting agents in a society 
(or among societies) [^[|. Culture is here defined as a set of 
individual attributes that are subject to social influence. The 
agents are placed in the nodes of a square lattice, and the indi- 
vidual attributes of each site are defined by a set of F features, 
each taking one of q possible traits. The dynamics of social 
influence takes into account that "the more similar an actor is 
to a neighbor, the more likely the actor will adopt one of the 
neighbor's traits" [Bp. Specifically, two neighboring sites cho- 
sen at random interact with a probability proportional to the 
number of common features they already share. An interac- 
tion translates into switching one of the different features of 
one site to the same trait of the neighboring site. In princi- 
ple, this dynamics tends to homogenize neighbors, and hence 
to reach a completely ordered state (a monocultural state) in 
which each feature has the same trait in all lattice sites. What 
Axelrod emphasizes about his model is that it illustrates how 
the mechanism of local convergence can generate global po- 
larization: He finds that the system gets trapped in multicul- 
tural states with a number of different stable homogeneous do- 
mains, called cultural regions. The number of domains in this 
disordered state is used as a measure of cultural diversity. The 
quantitative Statistical Mechanics analysis partially modifies 
Axelrod's conclusion Defining an order parameter as the 
relative size of the largest homogeneous domain, one finds a 
nonequilibrium order-disorder phase transition controlled by 
the number of traits q, which measures the initial degree of 
disorder of a random configuration. Symmetry breaking lead- 



ing to dominance of a given culture occurs for q < q c , while 
the order parameter vanishes for q > q c . Therefore local con- 
vergence is efficient for q < q c and global polarization only 
persists for q > q c . 

The analysis of corresponds to a zero-temperature dy- 
namics in which fluctuations are neglected. Here we address 
the issue of the effects of random perturbations (noise) fol- 
lowing the original suggestion of Axelrod, "Perhaps the most 
interesting extension and, at the same time, the most diffi- 
cult one to analyze is cultural drift [modeled as spontaneous 
change in a trait]" [||Q. Axelrod also identified that a crucial 
question in this case is to determine a characteristic time scale 
of the dynamics. We find that this characteristic time scales 
with the system size. A ^-independent noise induced transi- 
tion [ [To| ] occurs when the noise rate r is of the order of the in- 
verse of this time scale T: for r < T~ x noise induces ordered 
monocultural configuration because the disordered multicul- 
tural states are metastable; for r > T -1 symmetry restoring 
by noise is efficient and the disordered multicultural state pre- 
vails. In the thermodynamic limit of large systems the condi- 
tion r > T^ 1 — > is satisfied for any finite r. In this limit, 
the fundamental idea of local convergence generating global 
polarization is recovered. 

The model we study is defined [BJ by considering N agents 
as the sites of a lattice. The state of agent i is a vector of F 
components (cultural features) (era, er^, • ■ ■ , <Jip). Each a if 
is one of the q integer values (cultural traits) 1, . . . , q, initially 
assigned independently and with equal probability 1/q. The 
time-discrete dynamics is defined as iterating the following 
steps: 

1. Select at random a pair of neighboring lattice sites 

2. Calculate the overlap (number of shared features) 

3. If < < F, the bond is said to be active and 
sites i and j interact with probability l(i,j)/F. In case 
of interaction, choose g randomly such that a lg ^ aj g 
and set a% g = Uj g . 

In any finite lattice the dynamics settles into an absorbing 
state, characterized by the absence of active bonds. Obviousy 
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all the completely homogeneous states are absorbing. Inho- 
mogeneous states consisting of two or more homogeneous do- 
mains interconnected by bonds with zero overlap are absorb- 
ing as well. 

First we analyze the stability of these absorbing states. Be- 
fore we turn to the original two-dimensional geometry of the 
model, we consider the one-dimensional lattice with nearest- 
neighbor interaction. For this geometry the total negative 
overlap V = — X^i* i + 1) is a Lyapunov function: a 
function that never increases during the dynamical process. 
Assume that at a given time step the interaction across an ac- 
tive bond (i, i + 1) switches the state of one of the features of 
site i. Then the overlap across that bond increases by one unit. 
However, the overlap across the other bond (i — 1, i) of site i 
can at worst decrease by one unit. Taking into account that all 
other terms in V do not vary, we find that in any interaction, 
V never increases. There is a multiplicity of q F ground states 
of V that correspond to the homogeneous (ordered) configu- 
rations a if = <jjf and /, where the potential V takes its 
absolute minimum V m i n = —NF. Any other (disordered) ab- 
sorbing configuration is not the absolute minimum of V and 
therefore is a metastable minimum of the Lyapunov potential. 

We have not found a parallel argument for e?-dimensional 
lattices (d > 1). However, by extensive numerical simulations 
we have verified the metastability of disordered configurations 
for d = 2 with nearest neighbor interaction as well. The ab- 
sorbing states are subject to single feature perturbations, de- 
fined as randomly choosing i £ {1, N}, / € {1, . . . , F} 
and s 6 {1, . . . , q} and setting a%f = s. Then the simulations 
are designed as follows: (I) Draw a random initial configu- 
ration. (II) Run the dynamics by iterating steps (1), (2) and 
(3), until an absorbing state is reached. (Ill) Perform a sin- 
gle feature perturbation of the absorbing state and resume at 
(II). We find that by these cycles of relaxation (step II) and 
perturbation (step III) the system is driven to complete order, 
where cr,/ = <?jf, Vi, j and /. Thus, as in the one-dimensional 
case, only the completely ordered configurations are stable, all 
other absorbing configurations are merely metastable. In Fig. 
[j] we have plotted a typical evolution of the order parameter 
considered in [BR, that is, the size of the largest homogeneous 
domain S max . We observe that, although eventually its value 
could decrease, S max typically increases in each cycle of per- 
turbation and relaxation JlTp . 

We next consider the effect of cultural drift modeled as 
a random perturbation acting continuously on the system 
(noise). This is implemented by including a fourth step in 
the iterated loop of the model: 

4. With probability r, perform a single feature perturba- 
tion. 

Note that now the cultural imitation among sites and the per- 
turbations are two processes acting on the state vectors on time 
scales separated by the factor r. In contrast to the previous 
scenario, the system is not necessarily in an absorbing config- 
uration when a perturbation occurs. 



Figure || shows the average order parameter (S'max) as a 
function of the number of traits q. In the case r — 0, we 
observe the expected transition from order to disorder as q is 
increased |]^,[T2|]. Strikingly, no such transition is observed for 
noise rate r ^ 0. In fact, (5 max ) is practically independent 
of q except for small values. This becomes clear by notic- 
ing that with probability 1/q a perturbation does not change 
the configuration. Therefore considering the effective noise 
rate r' — r(l — 1/q) the data collapse onto a single curve 
(S'max 0"')) ( see Fig- §)• This curve identifies, for a fixed 
size of the system, a continuous order-disorder transition con- 
trolled by the noise rate. In addition, it shows that the transi- 
tion has universal scaling properties with respect to the value 
of q. 

It is worth noting that the limit r — > does not recover the 
results for the original model without noise, corresponding to 
r = 0. The fact that for r^Owe find an ordered state for all 
values of q is linked to the metastable nature of the inhomo- 
geneous absorbing states previously discussed. Perturbations 
at a vanishingly small rate are sufficient to allow the system to 
escape from these states. 

The order-disorder transition is also observed in the dis- 
tribution of sizes of homogeneous domains. For small noise 
rates typically a single cluster spans the whole system. As the 
noise rate is increased smaller clusters become more and more 
abundant. At an intermediate value of the noise rate the dis- 
tribution of domain sizes follows a power law. The exponent 
is approximately —2. 

What causes the onset of disorder with increasing noise 
rate? For a sufficiently small noise rate the system has enough 
time to relax to an absorbing configuration between perturba- 
tions. Then the situation would be comparable to the previ- 
ously studied case of alternating perturbation and relaxation: 
after a transient of reducing the disorder from the initial con- 
dition like the one shown in Fig. 1, the system will spend 
most of the time in one of the homogeneous configurations. 
At variance with the case of alternating perturbation and re- 
laxation where the ordered state cannot be left once reached, 
there will be, on long time scales, random jumps taking the 
system from a homogeneous state to another homogeneous 
one. If the noise rate is increased such that the typical time 
1/r between perturbations is shorter than the average relax- 
ation time, perturbations are "accumulated" in the system and 
disorder is built up. According to this picture the average re- 
laxation time T of perturbations of a homogeneous state sets 
the transition where rT = Oil). 

An approximate argument for the calculation of T is as fol- 
lows. A single feature perturbation of an ordered state at time 
t = induces a "damage" of size x[t = 0) = 1 in one of 
the components. In the following time steps the damage may 
spread until an ordered state is reached again by x(t) = or 
x(t) = N. The probability D x for x to increase or decrease 
is the fraction of active bonds, which again depends on x but 
also on the current shape of the damage cluster. This com- 
plication is avoided in a mean-field description, where a bond 
exists between any two sites. Then the number of active bonds 
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is simply the number of pairs of sites i carrying the dam- 
age and sites j not carrying the damage. Given damage size x 
we find x(N — x) active out of N(N — l)/2 bonds, yielding 
an update probability (for N > 1): D x = 2N~ 2 x(N - x). 
The average first exit time, that is the average number of time 
steps t x required to reach an absorbing boundary (x = or 
x = N) when starting from x(t = 0) = x, fulfills the recur- 
sion relation [ 13 1 



t x = ~D x (t x ^i + T x+ l) + (1 - D x )t x + 1 . (1) 

with the boundary conditions tq = tn = 0. This equation has 
the solution 
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The average relaxation time T is the average first exit time 
for an initial damage size x = 1. Thus we obtain T = 
n = N (7 + ln(iV - 1) + 0(7V -1 )), with Euler's constant 
7 w 0.577. For large iV we approximate 



T = iVln(jV), 
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so that the average relaxation time of perturbations diverges 
with increasing system size. 

This result, together with our argument that the transition 
occurs for rT = 0(1), suggests to consider systems of dif- 
ferent size. We always find a smooth transition from the or- 
dered to the disordered state as r is increased, see Fig. || The 
noise rate necessary to induce disorder is seen to decrease with 
growing system size. Fig. ^ confirms that T is the relevant 
time scale for the onset of disorder: The data of Fig. Q collapse 
into a single curve when plotted as a function of a rescaled 
noise rate R = rT = rN ln(iV) which incorporates noise rate 
r and system size N. Recalling our result for a fixed system 
size, we conclude that there is a universal scaling form for the 
order parameter describing the transition: (S ma , x (r, q, N)) = 
(S max (fir,q,N)) with f(r,q,N) = r(l - l/q)N\n(N). 

We finally note that for the correct prediction of the charac- 
teristic time scale T and the change of behavior with noise rate 
r, the dependence of the interaction probability on the sites' 
overlap need not be taken into account. This motivates the in- 
troduction of a "decoupled" version of the dynamics in which 
a site always adopts the trait of the chosen neighboring site in- 
dependently of the number of shared features. It turns out that 
the size of the largest cluster follows the same dependence on 
the parameters q, r > and N as in the original model (inset 
of Fig |]). This shows that in the presence of noise our main re- 
sults are insensitive to one of the basic premises of Axelrod's 
model. 

In summary, we have described the scaling properties of an 
order-disorder transition induced by noise which occurs at a 
size-dependent value of the noise rate T^ 1 (N). For a finite 
system and r < T^ 1 ordered monocultural configurations are 
induced by noise because disordered multicultural configura- 
tions are metastable. In the thermodynamic limit TV — > 00, 



the condition r > T^ 1 (N) — > is satisfied for arbitrarily 
small noise rate and the system remains in a disordered or 
multicultural state. In any case, cultural drift changes the na- 
ture of the ordered and disordered states. In the ordered state 
the system is not stuck in a single homogeneous configura- 
tion, but during long time scales visits in succession a series 
of monocultural configurations. Likewise, the disordered sate 
is not a frozen configuration but an evolving state with noise 
sustained dynamics. Thus, the cultural drift appears to be a 
relevant variable which drastically modifies the dynamics of 
Axelrod's cultural model. 
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FIG. 1. Metastability of disordered absorbing configurations: 
The order paramater S max normalized with systemsize N is shown 
as a function of perturbation cycles for a square lattice of N = 100 2 
sites using F = 10 and q = 60. 
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FIG. 2. (Smax) as a function of q for different values of r. The 
dashed line is for r — with a transition at q = q c « 50. The simu- 
lations used AT = 50 2 sites with F = 10 features. Values (S max ) in 
Figures ^.^Q. and ^| are averages over 500 configurations. Measure- 
ments were taken after a relaxation time of 100 iV 2 time steps. 
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FIG. 4. (S ma x) as a function of noise rate for different values of 
N. Simulations were performed with F = 10 and q = 100. 
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rescaled noise rate R = r N ln(N) 
FIG. 5. (Smax) plotted as a function of the rescaled noise rate 
R = rNlnN. Displayed symbols and model parameters are the 
same as for Fig. ^ The inset is the corresponding plot for the decou- 
pled model, see the main text. 
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effective noise rate r' 
FIG. 3. (Smax) as a function of the effective noise rate 
r' — r(l — 1/q) for different values of q. Simulations have been 
run in systems of size N = 50 with F = 10. 
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